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Lepton flavor violating Higgs boson decays (LFVHD) are studied in the context of the Minimal 
OO , Supersymmetric Standard Model (MSSM) being embeded in a mSUGRA scenario that is 

£NJ ■ enlarged with three right handed neutrinos and their supersymmetric partners, and with the 

neutrino masses being generated by the seesaw mechanism. We compute the partial widths 
for these decays to one-loop order and analyze numerically the corresponding branching ratios 
in terms of the mSUGRA and seesaw parameters. We analyze in parallel the lepton flavor 
changing lj — > k-f decays and explore the maximum predicted rates for LFVHD, mainly for 
H°,A° — > Tp, decays, by requiring compatibility with neutrino and BR(lj — » lij) data. We 
1 find LFVHD ratios of up to 10 -5 in some regions of the MSSM-seesaw parameter space. 
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1 Introduction 



The observed neutrino masses do require a theoretical framework beyond the Standard Model 
of Particle Physics with just three massless left-handed neutrinos. Within the MSSM-seesaw 
context, which will be adopted here, the MSSM particle content is enlarged by three right 



handed neutrinos plus their corresponding supersymmetric (SUSY) partners, and the neutrino 
masses are generated by the seesaw mechanism. Three of the six resulting Majorana neutrinos 
have light masses, m Ui ,i = 1, 2, 3, and the other three have heavy masses, m^, i = 1, 2, 3. These 
physical masses are related to the Dirac mass matrix mp, the right-handed neutrino mass matrix 
rriM, and the unitary matrix Umns by diag(m Ul ,m U2 ,m U3 ) ~ C^ Arg (— mom^- m^)?7MJVS and 
diag(mN 1 ,mj^ 2 ,m]y 3 ) ~ rriM, respectively. Here we have chosen an electroweak eigenstate basis 
where tum and the charged lepton mass matrix are flavor diagonal, and we have assumed that 
all elements in wid = Y v < H2 >, where Y u is the neutrino Yukawa coupling matrix and 
< H2 >= v sin/? (v = 174 GeV), are much smaller than those of mu- The two previous 
relations can be rewritten together in a more convenient form for the work presented here as, 
mj) = im c ^ a91 ^ 2 Rmi ta9l ^ 2 U^ NS , where R is a general complex and orthogonal 3x3 matrix, 
which will be parameterized by three complex angles 6i, i = 1,2,3. 

One of the most interesting features of the MSSM-seesaw model is the associated rich phe- 
nomenology due to the occurrence of lepton flavor violating (LFV) processes. Whereas in the 
standard (non-SUSY) seesaw models the ratios of LFV processes are small due to the small- 
ness of the light neutrino masses, in the SUSY-seesaw models these can be large due to an 
important additional source of lepton flavor mixing in the soft-SUSY-breaking terms. Even 
in the scenarios with universal soft-SUSY-breaking parameters at the large energy scale asso- 



dated to the SUSY breaking Mx, the running from this scale down to wim induces, via the 
neutrino Yukawa couplings, large lepton flavor mixing in the slepton soft masses, and provides 
the so-called slepton- lepton misalignment, which in turn generates non-diagonal lepton flavor 
interactions. These interactions can induce sizable ratios in several LFV processes with SM 
charged leptons in the external legs, which are actually being tested experimentally with high 
precision and therefore provide a very interesting window to look for indirect SUSY signals. 
They can also induce important contributions to other LFV processes that could be meassured 
in the next generation colliders, as it is the case of the MSSM Higgs boson decays into r/2, re 
and jj,e which are the subject of our interest. 

Here we compute the partial widths for these lepton flavor violating Higgs boson decays 
(LFVHD) to one-loop order and analyze numerically the corresponding branching ratios in 
terms of the mSUGRA and seesaw parameters, namely, Mq, Mi/%, tan/3, and R. For the 
one loop running of the parameters we use the mSUSPECT programme. We analyze in parallel 
the lepton flavor changing lj — > ^7 (i ^ j) decays and explore the maximum predicted rates for 
LFVHD, mainly for H°,A° — > r/2 decays, by requiring compatibility with BR(lj — > ^7) data. 
For these we use the present experimental upper bounds given by^ \BR(/j, — > ej)\ < 1.2 x 10 -11 , 
\BR(t -> < 3.1 x 1(T 7 and \BR(t -> &y)\ < 2.7 x 10~ 6 . 

For the numerical analysis we choose, Mx = 2 x 10 16 GeV and Aq = 0. The Umns 
matrix elements and the m Ui are fixed to the most favored values by neutrino data ^ with 

J&r^i = 0.008 eV, ^Am 2 atm = 0.05 eV, 9 12 = 9 sol = 30°, 9 23 = 9 atm = 45°, 9 13 = 0° and 
5 = a = [3 = 0. We consider two plaussible scenarios, one with quasi-degenerate light and 
degenerate heavy neutrinos and with m Ul =0.2 eV , m U2 = m ui + 2m so1 , m V3 = m vi H — 2m atm 
and TTT-TVi = rriN 2 = mjsf 3 = misf, and the other one with hierarchical light and hierarchical heavy 
neutrinos, and with m ui ~ OeV ,m U2 = y Am^ , m Vi = yj Am^ tm and < < itin 3 - 

This is a reduced version of our more complete work^ to which we address the reader for 
more details. 



2 Numerical results and conclusions 



We show in figs. Q through Q the numerical results for the branching ratios of the LFVHD 
together with the branching ratios for the relevant lj — > ^7 decays. The results of BR{Hq — > r/i) 
as a function of m^, for degenerate heavy neutrinos and real R, are illustrated in fig. (J^i), for 
several tan/? values, tan/3 = 3, 10, 30, 50. Notice that in this case, the rates do not depend on R. 
The explored range in rajy is from 10 8 GeV up to 10 GeV which is favorable for baryogenesis. 
We also show in this figure, the corresponding predicted rates for the most relevant lepton decay, 
which in this case is fi — > cy, and include its upper experimental bound. We have checked that 
the other lepton decay channels are well within their experimental allowed range. The ratios 
for Aq decays, not shown here for brevity, are very similar to those for Hq decays in all the 
studied scenarios in this work. We have also found that the ratios for the light Higgs boson, ho, 
behave very similarly with m^v and tan (3 but arc smaller than the heavy Higgs ones in about two 
orders of magnitude. From our results we learn about the high sensitivity to tan (3 of the LFVHD 
rates for all Higgs bosons which, at large tan (3, scale roughly as (tan (3) , in comparison with the 
lepton decay rates which scale as (tan (3) 2 . The dependence of both rates on mjv is that expected 
from the mass insertion approximation, where BR(H X — > Ijk), BR(lj — > ^7) oc \rriN log(mAr)| 2 . 
We find that the largest ratios, which are for #0 an d Aq, are in any case very small, at most 
10~ 10 in the region of high tan/3 and high m^. Besides, the rates for \i — > e7 decays are below 
the upper experimental bound for all explored tan/3 and values. The branching ratios for 
the Higgs boson decays into re and \xe are much smaller than the r/2 ones, as expected, and 
we do not show plots for them. For instance, for = 10 14 GeV, and tan (3 = 50 we find 
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Figure 1: Dependence of LFV ratios with mjv (GeV) for degenerate heavy neutrinos and real R and for several 
values of tan/3, (la) BR(Ho — > rp,). (lb) BR(n — > ej). The horizontal line is the experimental upper bound. 
In both plots, the solid, dashed, dashed-dotted and dotted lines are the preditions for tan/3 = 3, 10,30 and 50, 

respectively, and M = 400 GeV, M 1/2 = 300 GeV. 

BR(HW -» Tpt)/BR(HW Te) = 4 x 10 3 and BR{H&> -» rp)/BR(H^ -> p,e) = 1.2 x 10 6 
for the three Higgs bosons. 

The case of hierarchical neutrinos gives clearly larger LFV rates than the degenerate case, 
as can be seen in figs. (|2|>. © and (@J). However, we will get restrictions on the maximum 
allowed Higgs decay rates coming from the experimental lepton decay bounds. For instance, 
the case of real 9i, that is illustrated in fig. (j2J shows that compatibility with fx — > e^f data 
occurs only in the very narrow deeps at around 9\ = 0, 1.9 and tt. Notice that it is precisely at 
the points 9\ = 0,tt where the BR(Hq, Aq — > rp) rates reach their maximum values, although 
these are not large, just about 10 -8 . We have checked that for lower tan/? values, the allowed 
regions in 9\ widen and are placed at the same points, but the corresponding maximum values 
of the LFVHD rates get considerably reduced. For the alternative case, not shown here, of real 
9i 7^ 0, with 9\ = #3 = we get a similar behaviour of BR{H X — > rp) with 92 than with 9\, and 
the maximum values of about 10 -8 are now placed at 9i = 0, tt. In contrast, BR(fi — » erf) is 
constant with 9% and reach very small values, well below the experimental bound. In particular, 
for tan (3 = 50, M = 400 GeV and M 1/2 = 300 GeV it is 10" 19 . Regarding the dependence with 
#3, not shown here either, a reverse situation is found, where BR(H X — > rp) is approximately 
constant and, for the heavy Higgs bosons, it is around 10 -8 . On the contrary, BR(fi — ► ey) 
varies but it is always well below the experimental upper bound. In addition, we have checked 
that the BR{t — > ^7) and BR(t — > erf) rates are within the experimental allowed range in 
all cases. In conclusion, for real R we find that the maximum allowed LFVHD rates are at or 
below 10 -8 . The case of complex R is certainly more promissing. The examples illustrated in 
figs. (jHJ) and (J3J) are for the most favorable case, among the ones studied here, of complex #2 7^ 
with 9\ =$3 = and show that considerably larger BR(H X — > rp) rates than in the real R 
case are found. Regarding the dependence with 92, we find that for the explored values with 
(\92\,Arg(92)) < (3.5, 1), the Higgs rates grow with both |^| and Arg{92) and, for the selected 
values of the MSSM-seesaw parameters in fig. ©, they reach values up to around 5 x 10~ 5 . We 
have checked that the predicted rates for BR(fi — > e~f) are well below the experimental upper 
bound, being nearly constant with 92 and around 10 -19 . Similarly, for the r — > decay. Notice 
that the smallness of these two decays, in the case under study of 92 ^ 0, is not maintained if our 
hypothesis on #13 = is changed. For instance, for #13 = 5°, which is also allowed by neutrino 
data, we get BR(fi — > erf) ~ 1.8 x 10~ 8 well above the experimental upper bound. Therefore, 
in this case of complex 92 7^ 0, the relevant lepton decay is r — > ^7 which is illustrated in 
figs. |31) and together with its experimental bound. For the set of parameters chosen in 
fig. ©, we get that the allowed region by r — > fij data of the (\02\, Arg(92)) parameter space 
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Figure 2: LFV ratios for hierarchical heavy neutrinos and real R with Qi ^ 0, 82 = 63 = and (tun-i , mjv 2 , mjv 3 ) = 
(10 8 ,2 x 10 s , 10 14 ) GeV. Here, tan/3 = 50, M = 400 GeV and M 1/2 = 300 GeV. (2a) BR(H X -> 774). 
Solid, dashed and dashed-dotted lines (the two later undistinguishible here and in all plots) correspond to 
H x — (ho, Ho, Ao) respectively. (2b) BR()i — + e-y). The horizontal line is the experimental upper bound. 



implies a reduction in the Higgs rates, leading to a maximum allowed value of just 5 x 10 . 
The dependence of the LFV ratios with Mq and Mi/ 2 for hierarchical neutrinos are shown in 
fig. ((I}). We see clearly the different behaviour of the LFVHD and the lepton decays with these 
parameters, showing the first ones a milder dependence. This implies, that for large enough 
values of Mq or M±/ 2 01 both the BR(t — > ^7) rates get considerably suppresed, due to the 
decoupling of the heavy SUSY particles in the loops, and enter into the allowed region by data, 
whereas the BR(Hq — > 77Z) rates are not much reduced. In fact, we see in figs. and (|3f) 
that for the choice Mq = M\ji the r decay ratio crosses down the upper experimental bound 
at around Mq = 1100 GeV whereas the Higgs decay ratio is still quite large ~ 4 x 10 -6 in 
the high Mq region, around Mq ~ 2000 GeV. This behaviour is a clear indication that the 
heavy SUSY particles in the loops do not decouple in the LFVHD. Notice also that it can be 
reformulated as non-decoupling in the effective H^Tfi couplings and these in turn can induce 
large contributions to other LFV processes that are mediated by Higgs exchange as, for instance, 
r — > /i/i/z. However, we have checked that for the explored values in this work of Mq, Mi/%, 
tan j3, R and mjVj that lead to the anounced LFVHD ratios of about 4 x 10 -6 , the corresponding 
BR(t — > fifJ-n) rates are below the present experimental upper bound. 

In summary, after exploring the dependence of the LFVHD rates with all the involved MSSM- 
seesaw parameters, and by requiring compatibility with data of the correlated predictions for 
fi — ► cy, t — > e^/ and r — > fi'y decays, we find that BR{Hq, Aq — > rp) as large as 10~ 5 , for 
hierarchical neutrinos and large Msusy in the TeV range can be reached. These rates are close 
but still below the expected future experimental reach of about 10 -4 at the LHC and next 
generation linear colliders. 
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Figure 4: Dependence with M (GeV) and M 1/2 (GeV) for (m Nl ,m N2 ,m N3 ) = (10 s , 2 x 10 s , 10 14 ) GeV, <9i = 
6» 3 = and tan/3 = 50. (4a) BR(H -> r/Z) versus M (GeV) for M 1/2 = 300 GeV and 9 2 = 7re°- 4 \ (4b) Same 
as (4a) but for BR(r -> ^y). (4c) Bi?(# -> r/t) versus Mi /2 (GeV) for M = 400 GeV and 82 = 7re°' 8 \ (4d) 
Same as (4c) but for BR(t -> ^7). (4e) Si2(JJ -> TjLt) versus M = M 1/2 (GeV) for 6 2 = 7re°- 8 \ (4f) Same as 
(4e) but for BR(r — * /17). The horizontal lines are the upper experimental bound on BR(t — > fi'y). 
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